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The radiation loss of small amplitude oscillons (very long-living, spatially localized, time depen- 
dent solutions) in one dimensional scalar field theories is computed in the small-amplitude expansion 
analytically using matched asymptotic series expansions and Borel summation. The amplitude of 
the radiation is beyond all orders in perturbation theory and the method used has been developed 
by Segur and Kruskal in Phys. Rev. Lett. 58, 747 (1987). Our results are in good agreement with 
those of long time numerical simulations of oscillons. 



I. INTRODUCTION 



Time dependent spatially localized solutions in various field theories, which are long living in the sense of staying 
localized much longer than the light crossing time, have been already found in the seventies [l], 0, [H, H[ and ever since 
they are still attracting considerable interest [H, @, 0, II] • An obvious reason for interest is the unexpected longevity 
of these objects, all of which exhibit nearly periodic oscillations in time. The presence of at least one real massive 
scalar field seems to be necessary in order that such long living, spatially localized field oscillations - oscillons - could 
form. Since oscillons appear in the course of time evolution starting from rather generic initial data, this provides 
another reason to consider them of physical importance. Oscillons have been found to form in physical processes, e.g. 
as a result of vortex-antivortex annihilation in two dimensional Abelian Higgs model Q , domain wall collapse in 4ft 
theory [Icj . QCD phase transition [ill ], or during symmetry breaking in three dimensional Abelian Higgs model [12 ] . 
Therefore there is some reason to believe that oscillons (or configurations close to them) influence the dynamics, and 
they play a role in phase transitions, cosmology and the dynamics of extended objects (cosmic strings, domain walls, 
etc.), see e.g. [3 Q 13 • Importantly oscillons have been also found in the bosonic sector of the Standard Model 
p7l.ll8l.[l9j. There are attempts for including fermion fields in the study of oscillons [2(|. Oscillons resemble breathers 
of the one dimensional sine-Gordon (SG) theory, with the important difference that unlike true breathers they are 
continuously loosing energy by slow radiation. An oscillon just like a breather possesses a localized "core" , but differs 
significantly in its "radiative" reg ion outside of the core. Oscillons have been observed in various spatial dimensions, 
from d = 1 up to d = 6 [2l|, |22j |. there is, however, a marked difference between oscillons in d < 2 and in d > 2. 
In dimensions d = 1,2 oscillons can be well described by an adiabatic time evolution of breather-like configurations 
with an ever decreasing amplitude and with an increasing frequency tending towards a limit determined by the mass 
threshold p3l. l24l. [25j . In higher dimensions there are various types of oscillons, they exhibit instabilities and their 
behaviour is more complex [3, [1, [26| . Oscillons have been studied on a 1 + 1 dimensional expanding background in 

Refs. mmi. 

In this paper we consider oscillons in scalar field theories with a general self-interaction potential in one spatial 
dimension. So far most work on oscillons has been either purely numerical or based on various approximations. 
E.g. oscillon e nerg y and lifetime has been recently estimated in Ref. (29| . Our starting point is the small amplitude 
expansion [30L l31| which yields breather-like configurations with spatially localized cores. The small amplitude 
expansion yields an asymptotic scries for the core. Solutions of the field equations are either periodic in time or are 
radiating. In the first case standing wave tails are present outside of the breather core, whereas in the radiative case 
there are outgoing waves. These latter correspond to oscillons. Following Segur and Kruskal [32| and adapting the 



approach developed in Refs. [331 ] . 34j | we compute analytically one of the most important physical characteristics 
of oscillons, the amplitude of the outgoing wave responsible for their eventual demise. Although strictly speaking 
our results are only valid for small amplitude oscillons, this is not a major drawback, since in d = 1 the long time 
behaviour of any oscillon (which is what we are interested in) is determined by that of the small amplitude one. The 
energy loss of a small amplitude oscillon in. d = 1 to leading order can be written as 

£--*-". (i) 

where E is the energy of the oscillon core, A and B are constants determined by the theory. The corresponding 
equation for the </> 4 theory has been found first in Ref. [12] • Geicke [35[ has verified numerically the asymptotic energy 
loss following from Eq. (TTJ), i.e. E(t) ~ B/\nt. Since this is not a completely straightforward numerical exercise (one 
has to prepare very good initial data, one needs long time simulations to high precision, etc.) we have also made a 
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detailed numerical investigation. For small enough amplitudes, e, one has E oc e and it is clear that the energy loss is 
beyond all orders in s. In this work we compute the constants A, B analytically using matched asymptotic expansion 
and Borel summation techniques, as well as by an independent numerical method. Also we have been able to give 
convincing numerical evidence for the validity of the radiation law |T]) by preparing good initial data and performing 
long time simulations. We have also checked the validity of our analytical and numerical results on the example of 
the sine- Gordon theory where A is exactly zero. 

The plan of this paper is the following: First, in subsection lll Al we present the essential points of the small amplitude 
expansion in a one dimensional scalar field theory and relate it to the Fourier expansion of time-periodic solutions. 
Then we analytically extend the solution of the mode equations to the complex plane. In subsection III Bl we relate 
the radiative tail of oscillons to an exponentially small correction to the asymptotic expansion. In subsection III CI 
the radiation law, Eq. ([1]) is derived. In subsection III Dl we determine the amplitude of the exponential correction 
by solving numerically the Fourier mode equations in the complex plane. In subsection III El the Borel summation 
procedure is used to calculate this correction analytically. In second half of the paper numerical simulations are 
presented supporting the earlier analytical results. In subsections MI Al and IIIIBI we introduce the methods used in 
numerical simulations and for generating good initial data starting from the small amplitude expansion and employing 
a tuning respectively. In subsection IIII CI we estimate the lattice effects by putting the sine-Gordon breather on the 
lattice and measuring its energy loss. In subsection llll Dl wc investigate 4 oscillons evolving from initial data extracted 
from the small amplitude expansion. We found that the functional form of the theoretical radiation law, Eq. (JTJ) fits 
our data points; however, the fitted value of A differs from the theoretical one. One reason for the disagreement is the 
exceptionally small value of A in the 4 theory. Furthermore we have performed simulations starting with the initial 
data of Ref . Ref . (3f| . In subsection IIII El we investigate oscillon radiation in a specific 6 theory where the value A 
is maximal. In this case our numerical results agree with the theoretical one with satisfactory precision. 



II. ANALYTIC APPROACH TO OSCILLON RADIATION 

A. Small amplitude expansion 

We consider a real scalar theory in a 1 + 1 dimensional Minkowski space-time, with a general self-interaction 
potential, £/(</>). The equation of motion is just a non-linear wave equation (NLWE) given by 



>,tt 



= [/'(0H0 + ^.9 fc fc , (2) 



fc=2 



where <j> is a real scalar field. In Eq. ([2]) the mass of the field is chosen to be 1, and it has been assumed that the 
potential, U(<f>), can be written as a pow er series in <p where the are real constants. 

As explained in detail in Refs. [30l. l3ll| small amplitude solutions of Eq. © can be represented by the series: 

oo 

cb = J2^Mr,0, (3) 

k=l 

where e is a small parameter and the coordinates have been rescaled as r — ut and £ = ex, with lu = s/l — s 2 . 
The time dependence of the functions <j)k (r, Q is found to be determined recursively by a set of (forced) oscillator 
equations. For example the first few 4>k can be written as: 

01 = Pi (0 cos T 

02 = i 9 2P?(C)(cos(2r)-3) (4) 

03 = Pa(0 cost + ^_(4 5 2 - 3A)p?(C) cos(3t) , 

where Pi{Q, PsiCl are given in terms of the parameters of the potential and of a single function, 5(£) as: 

1 , (5) 



-A 2 - -\g\ + -9 5 - -g 294 + -g\\ Z(() - -A 5 2 5(C)(32 + 195 2 (C)) 



*® = ^> ?3(C) = ^ 

together with Z = 5(4 — 5 2 )/3, and A = 5g 2 /6 — Sg^/A > 0. The function 5(C) is a globally regular solution of the 
following nonlinear equation: 

d^5 



5 + 5 J = . (6) 
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By fixing the center of symmetry at £ = the regular solution of Eq. ([6]) is given simply by 



5(C) = %/2scch(C) 



(7) 



The only condition on the potential ensuring the existence of an exponentially decreasing solution to Eq.® is the 
positivity of the parameter A. It is clear that each term in the scries Eq. ([3]) is exponentially decaying in space and is 
periodic in time, hence if it would converge it would yield exponentially localized breathers. It has been demonstrated 
by Segur and Kruskal |32j for the example of the </> 4 theory with U{4>) = (j) 2 ((f> — 2) 2 /8, that the small amplitude 
expansion (J3j> does not converge, actually it is an asymptotic series. Nevertheless for sufficiently small values of e this 
asymptotic series constitutes an excellent approximation for oscillons of frequency u>(e) over a large interval in £ [3l| . 
The physical reason for the absence of spatially localized, exactly time-periodic breathers is simply that for generic 
potentials, time-dependent solutions of the NLWE @ radiate. This makes it very plausible that the asymptoticity 
of the series ([3]) in general models is due to radiative phenomena. An interesting prototype exception, i.e. when ([3]) 
converges, is the celebrated breather solution of the sine-Gordon (SG) theory, where U(<f>) = 1 — cos((p). The SG 
breather can be written as 



(x, t)s = 4 arctan 



w(e) cosh(ea;) 



4 arctan 



s(0< 



(8) 



It is a simple matter now to show that the small amplitude expansion of the SG breather yields a convergent series, 
indeed. Let us note here that in Eq. |(5J) |e| < 1, and it must not be small. By performing the e expansion of (JSJ) we 
obtain 



,Qb =V2 sS(0 [2 + e 2 - £ 2 S 2 (C)/4] sinr + V2e 3 S 3 {() sin(3r)/6 + 0(e 5 ) 



(9) 



and it is easy to see that this series converges for \sS\ < 1, thus for e < 1/V2 the small amplitude expansion of the 
SG breather converges. 

For a given theory the series ([3]) is unique. Since all of its terms are time-periodic, it represents a family of breather- 
like solutions in the sense of an asymptotic series. At this point it is natural to look for time-periodic solutions of the 
NLWE Eq. by expanding the field (j> in Fourier series: 



0(C,r) = $> fc (0 cos(fcr), 

fc=0 



(10) 



leading to an infinite set of mode equations for the <j>^ 
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6m,-p-q- A remarkable simplification takes place, if the potential is 



where 6 m ,± p ± q = 6 m ,p+q 

symmetric around zero, i.e. g2k = for k = 1, 2, . . . . In this case only the odd Fourier coefficients are nonzero in the 
Fourier expansion and the mode equations take the form: 



dC 2 



/ 2 2 

in lo 



4>, 



93 

4 



oo 

E 

m,p,g— 1 



5 n> ± m ±p±q 



n, m,p,q = 1, 3, 5 . 



(12) 



Equations (|11| admit solutions with a spatially well localized core and an oscillating (standing wave) tail whose 
amplitude tends to a constant for |x| — > oo. The asymptotic tail can be approximated as consisting of a superposition 
of standing waves of frequencies ncu, n = 2, 3, . . . . We note that for bounded solutions all modes $„ for n > 1 contain 
two parameters which can be interpreted as an amplitude and a phase of the corresponding frequency standing wave. 
We are interested in solutions for which the amplitudes of the tails are much smaller than that of the core, and these 
are the ones related to oscillons. Because of the existence of the asymptotic standing wave tail these solutions are not 
unique. Intuitively it is clear that for a given frequency, solutions with the smallest possible amplitude standing wave 
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tail should be close to the "inner part" of oscillons. It is quite plausible that for a fixed value of e, the solution with a 
minimal amplitude tail and being symmetric with respect to the origin, is actually unique. This hypothesis is supported 
by the results of Ref. j2^|, where such breather-type solutions have been named quasi- breathers. Another type of 
quasi-breather is the (most likely unique) solution, for which all modes, $„ — * exponentially for say x — > — oo, with 
a small amplitude oscillating tail in the other direction. We shall refer to such objects as asymmetric quasi-breathcrs. 

For sufficiently small values of e the tail amplitudes become exponentially small in e. In this case the core can be 
treated separately from the tail, and a linear superposition of the two is a very good approximation to the solution. 
One can verify that the small amplitude expansion of pip reproduces the terms of the asymptotic expansion (U). We 
remark that the small amplitude expansion represents an essentially unique core. As already mentioned the amplitude 
of the standing wave tail is C(exp(— 1/e)), i.e. it is beyond all orders in perturbation theory. Segur and Kruskal (SK) 
[32| has worked out a method to compute this "transcendentally small" amplitude on the example of the </> 4 theory 
using matched asymptotic series expansion. In the following we shall use the SK method to find the amplitude of the 
standing wave tail of quasi-breathers when e — > 0. The main idea of the SK method is to define an "inner" problem 
in the complex £-plane in the neighborhood of the singularity of S closest to the real axis. Clearly it is located at 
( = ±iir/2 where the function S has a simple pole. In fact close to the pole at in/ 2 

S(y) = ^ + ^ + 0((eyf) , (13) 
ey o 

where the rescaled variable y is defined as 

C = in/2 + ey. (14) 

The inner problem is defined by the rescaled variables near the singularity and keeping only the leading terms in e. 
For example, for symmetric potentials the inner equations are: 

d 2 . +- 



dy 2 



m,p,g— 1 



$n = J $ ™*P $ <3 S n, ±m±p±q + ■■■ (15) 



where it has been also used that uj 2 = 1 to leading order in e. We look for such solutions of the mode equations of the 
inner problem which can be matched to the solution of the outer problem. In our case the outer problem is defined by 
the analytic continuation of the small amplitude expansion from the real axis. The matching region is parametrized 
in the following way: 

{|ej/|«l (ey^O), |y|»l (|y| ->«>), -tt < arg(y) < --} . (16) 

Making use of the fact that the leading e order of <&„ on the real axis is proportional to e" S n it follows that in the 
matching region 

+00 (0) 

$ ° = Etf + <^ 2 ) ( 17 ) 

fe=2 y 
+00 (n) 

k=n y 

(n) 

where the coefficients a k are uniquely determined by the small amplitude expansion on the real axis. The first few 
terms in the expansion p7[)(|18[) . up to order 1/y 4 , can be written as 
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B. Correction beyond all orders 



In this subsection we will construct an exponential correction to the asymptotic series, Eq. (|19|) which after matching 
to the outer region determines the amplitude of the standing wave tail. A correction beyond all orders to a divergent 
series might seem meaningless at first glance, however, following Ref. [3^, we can give meaning to the method by 
finding a place where at least the imaginary part of the original series converges. This region is the imaginary axis 
Key = 0, where the algebraic asymptotic series (fill]) is real, so Im$„ = 0, and the imaginary part of the series 
converges trivially. We divide $„ into real and imaginary parts, 



(20) 



and decompose the inner version of Eq. (with e = and using the variable y) into real and imaginary parts. 

Then one can linearize the imaginary parts of the mode equations, and obtain coupled linear differential equations 
for f2 n along the imaginary axis. These equations contain fl m terms multiplied by various powers of ty^- As a first 
approximation, in the matching region these can be neglected, and one gets decoupled homogeneous linear differential 
equations with constant coefficients for Q n . For Qq the solutions are CIq = cxp(±y), which are oscillating, non- 
decreasing functions in the direction of the imaginary axis and have to be omitted, as they cannot be matched to the 
e expansion on the real axis. The solutions for fl\ are linear in y and the matching conditions forbid them as in the 
previous case. The solutions fl n for n > 1 which are tending to zero as Imj/ — > — oo are 



v n exp 



\-i\fi 



(21) 



Because of linearization one cannot determine the amplitudes v n ; methods to calculate them will be presented in the 
next subsections. 

If V2 7^ the dominant among the exponential corrections is O2 which will yield the leading term in the radiation. It 
is possible to get correction terms to v n by taking into account terms proportional to in the differential equations 
and substituting them by the leading order terms of the asymptotic expansion. Also considering the coupling between 
different f2„, we get 



Oo 



v 2 exp 
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2iV3(g% - 3A) 



cm 4, 



O I - exp 
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In the case of symmetric potentials <I > 2fcS are absent, so the dominant contribution comes from 



v 3 exp 
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exp 



(-i>/24y) 



(22) 



(23) 



Similar exponential correction appears in the neighborhood of the singularity —in/2. We should still match the 
correction to the imaginary part, Eq. (|22[) to the solution on the real axis. Hence we linearize the equation of <&2, 
(fTTj) about the quasi-breather core and get the following solution: 



S$ 2 = C sin 



in (\/3z 



(24) 



where C and a are arbitrary constants. We analytically continue <5<I>2 to the complex plane. We match it to the 
exponential correction (|22p obtained in the inner region around the pole wr/(2e), and to the corresponding expression 
around — iir/(2e). This determines C and a 



C = 2v2 exp 



V3tt 



a = . 



(25) 



Neglecting contributions from higher Fourier modes the derivative of the field <j> in the origin oscillates as 

73V 



d x (j)\x=o = 2v 2 exp -- 



2c 



v / 3cos(2f) 



(26) 



For symmetric potentials the above matching procedure works in a completely analogous way. The original problem 
of determining a periodic solution of the field equation, Eq. ([2]) is well posed if we impose boundary conditions. We 
use the boundary conditions introduced by Scgur and Kruskal in Ref. (32| . namely we require the field to vanish 
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at x — > — oo and that the solution remains bounded. These requirements provide sufficient conditions to make the 
solution unique. We obtained a time-periodic solution of the field equation which is asymmetric with respect to x = 0. 
ft is clearly not a breather, since it has a standing wave tail in x — > +00 and in this sense it is weakly localized. 

We can add a transcendentally small standing wave to $2 as it solves the corresponding linearized mode equation 
about the asymmetric quasi- breather (AQB) denoted by 4>aqb and it will give the leading order transcendental term 
in the new solution. Therefore one can write a symmetric configuration <f>s in the following form: 

4>S = 4>AQB + 4>st , (27) 

where <f> s t denotes a transcendentally small standing wave. <f> s t has to make <j)s symmetric about the origin. This 
requirement is enough to uniquely determine it. The derivative of 4>aqb at x = is given by (|26[) . hence 

/ V3tt\ cos (V3 x + £) / N 

<fs st = 2^ 2 exp - ^— ■ , cos (2t) , 28 

V 2e I sin £ 

where £ is an arbitrary constant. The symmetric configuration with the minimal amplitude tail, corresponding to 
£ = 7r/2, has been named quasi-breather in Ref. (3l| . This result is consistent with our knowledge about symmetric 
solutions, as when we take into account the standing wave tail in one mode, one of the two free parameters is required 
for symmetrization. From Eq. (|28|) one immediately gets the standing wave tail of the asymmetric quasi-breather: 

4>aqb = 4^2 cxp ^— ^ ^ sin ^\/3 x^j cos (2t) for x — > +00. (29) 

In conclusion we were able to determine the asymptotic field of the asymmetric and symmetric quasi-breathers up to 
one parameter (^3). In subsections III D l and III El we will determine this parameter. 

C. Radiation law for small amplitude oscillons 

We can construct a symmetric time dependent solution of the field equation by repeating the same steps as in the 
time periodic case and replacing the standing wave (j) s t with a moving wave 4> ra d- Let us denote the oscillon field with 
4> osc . Then from 4> osc = 4>aqb + <Ara<2 it follows that 

4>rad = ~2V2 CXp ( — ^ I sin ( V 3 X + 2t ) . 



2s 



- sin (y/3 x 



We find that the amplitude of the tail of the quasi-breather determined by £ = ir/2 is equal to the amplitude of the 
outgoing radiation from the corresponding oscillon. This property has been already noted in Ref. [26j . 

After the little digression on quasi-breather tails and the determination of the oscillon radiation field we focus on 
the radiation law. The asymptotic oscillon field from Eq. (f30|) is: 




2v n exp I sin ( y n 2 — 1 x — nt ) for x — > +00, (31) 



^ = - 2 „„e x p(-^i) si „(^, + n f ) fc^-co, ,32, 

where n = 2 for asymmetric potentials and n = 3 for symmetric potentials. 

The energy carried away by these oscillating tails determines the time-averaged radiation power W of the oscillon: 



^ = W = -AnV^lvl exp i^^^j ■ (33) 

Since the radiation field is transcendentally small, it is reasonable to assume that during its time-evolution the core 
of the oscillon goes through undistortcd quasi-breather states. This statement will be referred to as the adiabatic 
hypothesis. The energy content E of the core of the quasi-breathers as a function of e (or equivalently as a function 
of its frequency us) is easily determined, it is given by: 

E=^ + 0(e 3 ). (34) 
A 
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Now the equation determining the change of the core energy with time (energy loss) for oscillons can be seen to be 
of the form given by Eq. (TTJ) . From Eq. ((T|) one can easily deduce that the leading order late time behaviour of the 
energy is given as 

B ( 21nlnA . , 

In this subsection the radiation law for small amplitude oscillons has been determined up to a single unknown 
parameter, A. The problem of finding A or what is equivalent the parameters resp. will be done in the following 
two subsections. 



D. Determining the radiation amplitude by solving the complex mode equations numerically 

In this subsection we numerically determine the leading radiation amplitude coefficients, namely v-i for asymmetric 
potentials and v% in case of symmetric potentials. First, we consider the <j> theory, in which case the only non- vanishing 
coefficients in the expansion of the potential are g% = — § and g$ = | . We consider various order truncations of the 
Fourier mode equations pip in the region close to the singularity. To illustrate our method we present in more details 
the calculations for the simplest truncated system that radiates, i.e. that for which only up to cos(2r) modes are 
kept. For the inner problem, in the e — > limit, the mode equations are 



dy 2 



-*o = W - Wo + 7 (*? + ^2) (So - 1) + W^ 



dy 2 



+ -$1 ($ 2 - 2) (2$ + $a) 
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dy 2 



3 3 3 

3$ 2 = ^$0*2 + ^$1 ($0 + $2 - 1) - 3$ $2 + g$2 



Expanding the mode equations into powers of l/y 



(36) 



£<4" 



(37) 



we find that all coefficients are fixed after choosing a 2 and the sign in 
j < 2 and = for j < i. In order to agree as well as possible with (fT 



±2z/v3. Furthermore 
obtained from the e expansion evaluated 



af ) = 0, for 



at the singularity region we set a 2 = and a\ = —2i/y/3. Then it turns out that 



i+l+2j 



for all integer j. 



Raising the truncation order high enough by adding more Fourier modes to the system all coefficients will necessarily 
agree with (fi"9]) . The expansion consistent with (|36|) is 
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(38) 



According to expansion (|38[) the imaginary part CI2 = Im$2 vanishes to all orders on the imaginary axis. However, 
this is not a convergent but an asymptotic series, and the actual solution of (|36p may include an exponentially small 
correction to on the imaginary axis, in accordance with (j22|) . Of course, the value of ^ will depend on the chosen 
boundary conditions. The method introduced by Kruskal and Segur [32| is to integrate the differential equations 
numerically along a constant Imy = j/j line from a large negative Rey = y r value to the axis Rcy = 0. The boundary 
conditions at y = y r + iyi are given by the expansion (|38[) and its derivative, truncated to an appropriate order in 
l/y. This works well for modes with i > 1, but $0 has the tendency to exponentially blow up along constant Imy 
lines. This numerically problematic issue can be avoided by treating the equation for $0 as a two point boundary 
value problem, setting $0 — at the axis point y — iyi and using the asymptotic series as the boundary value at 
y = y r + iyi. The other modes are treated as initial value problems by specifying their value and first derivatives 
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at y = y r + iyi. We note that in case of symmetric potentials <f>2; = everywhere, and the integration procedure 
simplifies to pure initial value problem. In that case the first radiating mode is $ 3 . 

For the actual numerical integration of the (f> 4 system we have chosen various yi values in the interval [—5, —13], 
y r in [—20, —500], and the expansion in the initial data was truncated to orders from 4 to 20. The equations were 
generated by Maple and its default boundary value problem differential equation solver was used for the numerical 
integration. For y r < —50 the obtained values for were only changing in the fourth digits on varying y r . The choice 
of truncation order in 1 /y in the initial data was not changing the results in their less than fourth digits if the order was 
chosen larger than 8. In Tabled] we present the obtained yi dependence of the imaginary part f^, when y r = —300 
and the initial data is of order 15. According to (|22p we approximate to leading order as v£ = 2 exp (-\/3j/i) ■ 
The 1/y correction gives a more precise result z/g = ^2 cxp (\fiyi) /(l — l/(yiV3)). We note that for the special case 
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0.00935565 


0.00853442 


0.00863027 
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0.00911169 


0.00849838 


0.00853980 
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0.00895494 


0.00846614 


0.00848769 


0.00848446 
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8.31980 ■ 
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0.00884615 


0.00844008 


0.00845283 


0.00845151 



TABLE I: Dependence of the exponential correction Q.2 on the location of the integration line lray = yi for the cf> 4 theory 
truncated at cos(2r). The corresponding radiation amplitudes are a l so given by taking into account n-th order corrections 
in 1/y. 

of the 4> 4 theory the coefficient of the 1/y 2 correction in VI2 vanishes, and hence is also valid to second order. For 
larger values of \yi\ it is possible to improve the precision by adding even higher order corrections. 

From TablelH we may give a first estimate on the actual value of the radiation amplitude as = (8.45 ±0.03) • 10~ 3 . 
However this value turns out to change drastically when adding higher Fourier modes to the system (f3"6"l) . In Table 
HT1 we give the calculated values for v-i when keeping Fourier modes up to order cos(nr). It is somewhat surprising 



;;. 


V2 


2 


8.45 • 10 -3 


3 


-7.115 ■ 10~ :i 


4 


-8.431 ■ 10~ 3 


5 


-8.454 ■ 10~ :i 


6 


-8.454 • 10~ :i 



TABLE II: Dependence of the radiation amplitude on the order of the used Fourier expansions in the cj> 4 theory. 

that the addition of the cos(3r) mode changes the sign of v 2l while the magnitude is quite close to the proper value. 
The addition of higher than $6 modes does not make any significant change in the value of v-i- As a conclusion, we 
can state that v-i = (— 8.454 ± 0.01) • 10 -3 . The numerical value obtained by Kruskal and Segur differs by a factor of 
2 due to their use of complex notations for the Fourier modes. In our units their result is (—9.0 ± 2.0) • 10~ 3 . 

In case of symmetric potentials $i = for even i, and the numerical integration method simplifies considerably A 
specific example we have looked at is a specific (j) 6 theory, with C/'(0) = (f>— (p 3 + (f> 5 . When using the mode equations 
only for $1 and $3 we get V3 = —0.91026. Adding the fifth and seventh modes changes the result to —0.90982 and 
—0.90977, respectively. Since the addition of even higher modes do not make significant change, we can state that for 
the (f> 6 theory: 

v z = -0.9098 ±0.0001. (39) 

An important check of the reliability of our method is to calculate 1/3 for the sine-Gordon potential U(4>) = 1 — cos(0). 
In this case, when taking into account all the mode equations and all expansion coefficients in the expansion of 
the potential, the exact result is known to be zero. Instead of changing these independently, when we solved mode 
equations up to order $^ we assumed that gj = for j > i. As it can be seen from Table [TtTI . the results appear to 
tend to zero fast as i increases. 
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i 


U3 


4 


2.32 


6 


-0.2316 


8 


8.35 ■ 10" 3 


10 


1.14 ■ 10~ 5 


12 


7.2 • 10" 7 



TABLE III: Radiation amplitude v% for the sine-Gordon theory truncated to order i in both the mode equations and potential 
expansion. 



E. Determining the radiation amplitude by Borel summation of the algebraic asymptotic series 

In this subsection we will solve Eq. lfl5|) using the algebraic asymptotic series ansatz (fl8|) in the neighborhood of the 
singularity. Our considerations will be applicable to the case of symmetric potentials. At the end of this subsection 
we will briefly discuss the problem of asymmetric potentials. The solution is unique as we have to match it to the 
original asymmetric quasi-breather continued analytically from the real axis. We truncate both the Taylor expansion 
of the potential and the Fourier expansion in order to have a finite set of equations and will work until we reach 
convergence as in subsection III Dl 

We will demonstrate how the determination of coefficients works for big k in leading order of k in the case of 
minimal system $i and $3 with a cubic nonlincarity. Since and are vanishing for even k we redefine the 
coefficients in order to get a more convenient form: 

*i= l E^T (4°) 

k=l y 



*3 = *£-^T- (41) 



k=2 y 



We will show in the following that the behaviour consistent with Eq. ([15)) to leading order in k is 

2k 2 



~ (42) 



^-Ki-lf^^. (43) 

The constant K can be determined by solving the equations up to some large order k and matching the gained 
coefficients to the determined asymptotic behaviour. To do so, we write up the mode equations: 



i T $ 1 = ^S I ($; + 4 I 4 s + 2«|) (44) 



and then determine the equations for the coefficients of l/y 2k ~ 1 keeping only terms of order (2k — 4)!: 



$ 3 = j (*? + 6$?$ 3 + 3*1) , (45) 



{2k - 3)(2fc - 2) A k ^ + ^ AjB k ^ = (46) 
(2k - 3)(2fc - 2) B k ^ +8B k + ^- A\B k -i = . (47) 



It is easy to figure out the value of A\ from the matching conditions (fl"9|) . which gives A\ = —-J—8/ (3</3). (Obviously 
we get the same result with an indeterminate sign by solving Eq. (|4~4"|) for the coefficient of 1/y 3 -) Using the value 
of A\ we get the asymptotic behaviour of A k from Eq. (|4"6")) already given in Eq. (|4"2")) . From Eq. (|47p we can even 
determine the O (l/fe 2 ) corrections to the asymptotic behaviour of B k for large k as 

, k (2k-2)\ A 1 , 5 



^~^H)W [ 1 + k + W^ ■ (48) 
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This formula enables us to determine the numerical value of K very precisely from B^s with moderate k value. 

By examining the structure of the equations and making use of the fact that the algebraic asymptotic series of 
<f>„ starts only at l/y n it is easy to prove that the above asymptotics do not change. Even the O (-p-) corrections 
to -Bfc are not affected by the involvement of further modes or higher order nonlinearities. The only effect of the 
introduction of further Fourier modes and higher order nonlinearities which originate in the self-interaction potential 
is the changing of the value K . 

On the one hand, this result gives the proof of the asymptoticity of the series, i.e. the coefficients are not more 
gravely divergent than (2k — 2)!. On the other hand, this property allows us to Borel-sum the series. It will turn out 
that the behaviour of the Borel-summed series in the vicinity of its singularity gives us the dominant radiation field 
configuration, i.e. the exponentially small imaginary correction to $3 on the imaginary axis, in the matching region. 

The first step in the Borel summation is 



V(z) 



iB k 



JZk-X 



1)! 



OO 



^ 2k 

k=l 



2fc-l 



1 \VsJ 



-In 



1 + iz/V8 
l-iz/VS 



(49) 



This Borel summed series has logarithmic singularities at z 
Borel summed scries of $3 which we denote by $3 



±iy/8. The Laplace transform of V(z) will give us the 



$3(2/) 



dte"V 



(50) 



The integrand of Eq. ([50]) has logarithmic singularities at t/y = ±iy/8. It has been explained in Ref. 34 1 how to 
compute the integral. We only take into account the singularity t/y = i\f& and not t/y = —iV8, because we would like 
to get the correction in the imaginary part of $3 for y points with negative imaginary parts, as we aim to approach 
the real axis. For Imy < the other singularity stays away from the integration path, while the singularity from 
t/y = appears for t = iy/8y. In order to define the integral for Rey = we use the analyticity of $3. When y 
is in the matching region determined by Eq. (|16p the singularity in t is in the lower half-plane, i.e. the contour in 
Eq. (p)D|) is above the singularity. When Re y — > the singularity tends to the real axis and the contour must stay 
above the singularity. The logaritmic singularity of V (t/y) does not contribute to the integral and integrating on the 
branch cut starting from it yields the imaginary part 



Im$ 3 (y) 



_ t i Ktt i Kit 
dt e — - — = — - — exp 



(51) 



This result agrees with that of Eq. ([23]) . hence we were able to determine z/3 analytically: 



^3 



Ktt 



(52) 



We already gave the method for determining the value of K by solving linear equations recursively for the coefficients 
of the algebraic asymptotic series up to some large k values. In the first step, we will show that our method is consistent 
with the fact that the sine-Gordon breather does not radiate. We will truncate the SG potential in various orders 
and determine the value of K by solving the complex mode equations with the algebraic asymptotic series ansatz. 
We will experience a monotonous decrease in the value of K as we increase the order of the Taylor expansion. This 
can be interpreted in the following way. K determines the increase in the coefficients of the algebraic asymptotic 
series which correspond to the coefficients appearing in the small amplitude expansion. As we get closer to the SG 
theory the asymptotic series for small amplitude quasi-breathers in the truncated SG theory is less and less growing. 
On the other hand, this means that the corresponding oscillons would radiate slower, as the radiation amplitude is 
proportional to K. In the limit of SG theory we should get a convergent series in the small amplitude expansion 
and thus a non-radiating breather. Our numerical experiences show that we can even use a smaller number of mode 
equations than the biggest power in the Taylor expansion to reach satisfactory convergence for the value of K for the 
given theory. We collected the results for K in the truncated SG theory in Table IIVI 



In the second step will will be looking for a symmetric 6 theory in which oscillon radiation is the fastest, i.e. 
K takes the biggest value. We will use this theory in numerical simulations for the verification of the theoretical 
radiation law for oscillons because we hope to measure radiation rate in this theory the most accurately. We will 
make use of the fact that g 3 can be defined into the fields, thus the only essential parameter in a symmetric 4> e theory 
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Order of truncation 


K 


-^■num 


4 


1.486 


1.48 


6 


-0.1472 


-0.147 


8 


5.306 • 10 -3 


5.31 ■ 10" 3 


10 


7.306 • 10" 5 


7.2 ■ 10~ 5 


12 


4.686 • 10 -7 


4.6 • 10~ 7 



TABLE IV: The value of K in truncated SG theories from the solution of the same number of mode equation as the order of 
truncation. K denotes the value obtained by Borel summation, while -Knum is the result of the numerical solution of the mode 
equations on the complex plane. 



is 95- (.93 < in order to have a localized solution of Eq. ©.) To do so, we solved the theory for <J> 1 -$ 7 up to k = 30 
with the coefficient g§: 

m = l<f> 2 -\^+ 9 j ** (53) 

C/'(0) = 0-0 3 + 55 5 . (54) 

The value of K as a function of g$ can be found in Fig. [TJ We calculated K form B30 in the figure, however the 
exact value of k only matters in the sixth digit. We only deal with positive g$s as they are the ones with a stable 
vacuum. We see two zeros in this domain. This does not mean that we found breathers in the corresponding theories, 
these points represent oscillons the dominant radiation field of which is in fifth mode, $5. These configurations are 
extremely long living objects. The first zero is very close to the SG theory for which 175 = 3/10 with the 33 = —1 
normalization. Upon introducing further nonlincaritics this zero would move exactly to the point representing the 
SG theory. From Fig. [T| (75 = 1 seems to be a very comfortable choice for numerical simulation. We find the precise 
value of K for g$ = 1 to be: 

K = -0.57915 , K mm = -0.5792 , (55) 
where K num has been extracted from Eq. ([39]) . 




FIG. 1: The value of if as a function of gs for symmetric cj> theories 



In the end we will briefly discuss the case of asymmetric potentials the class to which the </> 4 belongs. In these 
theories the presence of the 4>o mode results in the following asymptotics for large ks in the <&o mode: 

«4 0) ~(2fc-2)!. (56) 

This behaviour dominates the asymptotics of all other modes as well. Because we do not have an alternating sign in 
this dominant behaviour, the Borel summed series will have a singularity on the real axis, hence we do not get an 
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imaginary correction to the asymptotic series from this calculation. Thus, we cannot determine the radiation this way. 
The asymptotics which would determine the radiation in the case of asymmetric theories is a? ~ (— l) fe (2k — 2)!/3 fe . 
This is not the leading behaviour of <f>2's coefficients and is significantly suppressed. We have not yet succeeded in 
determining the hidden alternating part of the coefficients, therefore we cannot determine the radiation amplitude by 
this analytic method. 



III. NUMERICAL SIMULATIONS 



A. Aspects of the numerical simulations 



The numerical simulation of oscillons and their decay were performed with a slightly modified version of the fourth 
order method of line code developed and used in Refs. [26|, Ha [13] • The spatial grid was chosen to be uniform in the 
compactified radial coordinate R defined by 

2R . . 

where K is a constant. By this choice the whole range < x < oo of the physical radial coordinate x is mapped to the 
interval < R < 1, avoiding the need for explicitly describing boundary conditions at some large but finite radius. 
In most of our simulations we worked with k = 0.05, which proved to be ideal for oscillons for which radiation was 
numerically observable. 

The frequency of oscillons is determined by measuring the time elapsed between two subsequent maximums of the 
field configurations in the origin. By integrating the energy density in an R = 0.873 (x = 146.7) sphere at every 
timeslice we determine the energy E of the oscillon core and trace the energy loss W = dE/dt. Although this way 
of calculating the energy contains an arbitrarity, for small amplitude oscillons the corresponding asymptotic field is 
exponentially small and the energy density of the radiation tail can be neglected compared to that of the core. We 
measure the radiation in a relatively short time interval from the slope of the energy as a function of time. 

Starting from initial data obtained from the small amplitude expansion we simulate the time evolution and radiation 
of oscillons. Although our analytic calculations are only valid for infinitesimal e, we find that they approximate oscillons 
quite well after a one-parameter tuning. 

Studying the slow energy loss of oscillons we determine a semi-empirical radiation law. This means that we keep 
the functional form (|33[) but fit the parameters appearing in it. We find satisfactory agreement with the theoretical 
value of the parameters obtained in the previous sections. We also follow the evolution of a single oscillon through a 
very long time interval. This process will be precisely described by the radiation obtained by the fit law proving the 
assumption made when making analytic considerations: the system evolves through undistorted quasi-breather states 
adiabatically. 



B. Initial data 



We used the asymptotic series ([3]) truncated to order N, at the moment of time-reflection symmetry as initial data 
for the numerical evolution code, 

N 

^=°) = £ £ ^=°>. (58) 
k=l 

The aim of our numerical analysis is to obtain oscillon states which arc as periodic as possible, which means that their 
basic oscillation frequency and amplitude changes very slowly, uniformly and monotonically. The time evolution of 
these "clean" oscillons can be approximated by adiabatic evolution through corresponding frequency quasi-brcathcr 
states. However, since the expansion ([3]) is not convergent, the initial data (|58|) differs from the intended quasi-breathcr 
configuration. In general, using it as initial data, first a small portion of the energy is quickly emitted by radiation, and 
then a very long living localized oscillating configuration remains. However, generally, the frequency and amplitude 
of this "unclean" oscillon state also possesses a lower frequency modulation. We observed that the amplitude of this 
modulation can be significantly decreased by multiplying the initial data by an appropriate constant. By making a 
one-parameter tuning code we were able to obtain oscillon states clean enough for studying their basic energy loss 
rate. Otherwise one could not distinguish between the energy emitted by the oscillon and the energy released by the 
decay of the low frequency modulation. 
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Through the whole domain of simulation the sum (|58j) with N — 3 proved to yield the cleanest oscillon states. It 
is worth mentioning, that even for quite large values of e, after the tuning the asymptotic series yields appropriate 
initial data, although the e value of the initial data and the one calculated from the frequency by e = yL — u> 2 during 
the time evolution may differ. 

We illustrate the main steps by the example of the c/> 4 theory. The first few terms of the asymptotic series are: 



6 (-o) 



(59) 



± I 125 4 103 2 
9 I ~ b + ~ b 



The most naive estimate for an ideal truncation of an asymptotic series is that we should find the order where the 
terms of the series are equally big and higher order terms are starting to grow from this threshold. From (|59p for 
0.3 < £ < 0.6 truncation of the series at third order appears to be appropriate. This is also supported by the results 
of the numerical simulations. In the (j> 4 theory we found that e = 0.6 oscillons are the biggest ones which we can 
clean from the noise and for which the adiabatic hypothesis works. Fig. [3] shows the obtained "unclean" oscillon 
states from various order initial data. It is apparent that the order N = 3 gives the state with the smallest amplitude 
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FIG. 2: Evolution from initial data in various orders for e = 0.6 in the cj> 4 theory; field value maximums at the center are 
plotted. 

modulation. As already noted, multiplying the initial data with a constant close to 1 decreases the amplitude of this 
modulation even more. Figure [3] shows how effectively this method smooths the oscillon for e = 0.5. 



C. Reliability of the numerical results concerning oscillon radiation 

The reliability of the subsequent results concerning oscillon radiation can be estimated by checking how well our 
code simulates the exactly periodic sine-Gordon breather. We use the breather field configuration ((5|) as initial data 
for the time evolution. We note that using the series (J9j) truncated at third order and employing the tuning method 
gives the same order of radiation. Table IVl contains for different resolutions the radiation power and W/ Eb where Eb 
is the total energy of the breather (Eb = 16e). The number of lattice sites are given by the formula Ni at = RES x 128. 
It can be observed that the duplication of the lattice points results in the decrease of radiation due to lattice effects 
with a factor of 30. From now on we will keep ourselves to the following rule: we only consider oscillons at a certain 
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RES 


W 




w 

E B 




4 
8 
16 
32 
64 


-1.18062 ■ 10 -6 
-3.66652 ■ 10" 8 
-1.14045 ■ 10 -9 
-3.5542 • 10~ n 
-1.10946 • 10" 12 


1.48 ■ 10~ 7 
4.58 • 10~ 9 
1.42 ■ 10 -10 
4.44 ■ 10~ 12 
1.38 ■ 10" 13 



TABLE V: Resolution dependence of radiation power W in the sine-Gordon theory. The third column gives the ratio compared 
to the energy of the breather Eb- The number of lattice sites are Ni at — RES x 128. 



resolution if the numerically calculated W/E is at least one magnitude bigger than the value given in Tabic [Vj If 
this condition fails, that means one has to use higher number of lattice sites to measure the radiation rate reliably, 
which, however, can make the simulation time impracticably long. This is the reason why the energy loss of oscillons 
for very small e cannot be obtained by our method. 



D. Oscillons in the </> 4 theory 



Confirmation of the adiabatic hypothesis and the radiation law 



In Table IVII we give the energy loss for the <f> 4 theory for the biggest possible interval where we are able to use our 
numerical approach. The lower bound is limited by the need of time to perform high resolution simulations; we can 
see only lattice effects below the lowest e values in the table. Above the upper bound the oscillon decay is very fast 
and we cannot determine W and e reliably. We will see the radiation law fail for large e, in these cases the interaction 
with the radiation field may become essential, and the system does not evolve through undistorted quasi-breather 
states. The V2 values in Table IVII are calculated from the measured radiation power and e by using the theoretical 
radiation law ([55]) . 

We intend to confirm the radiation law, Eq. (|33() numerically. We performed two fits with two free parameters on 
the logarithm of the data shown on Figure 0Ja). I n the first fit we took into account all the data points, while in the 
second we fitted for data points with e < 0.42 to get closer to the theoretical pole term. We define the semi-empirical 
radiation law: 



dE 



W 



-8V3 " 2 



v 2 • exp 



V^3 7T b 



(60) 



where v 2 and b are parameters to be fitted. Eq. (|60p is to be interpreted as follows. For finite values of e there are 
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TABLE VI: Radiation power W and e in the </> 4 theory from various initial data. For the oscillons e is measured from the 
frequency during time evolution. The column determines the resolution RES required for the simulation, while is 

calculated from Eq. (|33|l . 



various higher order corrections to the radiation law which are now encoded in two effective parameters, v 2 and b. 
These corrections originate in Eq. (f22|) . We remark that the value of btheory will increase accordingly. When fitting 
for all the data points the difference between b and b t h eory = 1 is bigger, than in the case of the fit for e < 0.42 points, 
because for smaller e values the correction due to other pole terms are smaller. 

If we set b = btheory in Eq. (j60|) we get back the theoretical radiation law, Eq. l|33|) . On Figure HJb) we plotted the 
v 2 values from Table IVTl We see that we are getting closer to the theoretical v 2 as e decreases. In Eqs. (|fJ5|) . we 
give the results of the fits of v 2 with b = btheory on all the data points and for data points with e < 0.42 respectively. 

The result of the first fit supports the functional dependence got from the theoretical arguments, while the second 
fit leads us closer to the theoretical values of the parameters: 

i/ 2 ju = 0.29 ± 0.03 bfu = 1.28 ± 0.2 for all data points, (61) 

v 2 jit = 0.12 ± 0.02 bfu = 1.17 ±0.2 for data points with e < 0.42, (62) 

v 2Jlt = 0.054 ± 0.004 b=l for all data points, (63) 

v 2Jlt = 0.034 ± 0.004 b = 1 for data points with e < 0.42, (64) 

V 2 ,theory = 0.00845 hheory = 1 ■ (65) 



Unfortunately the theoretical v 2 value of the theory is exceptionally small. This might be the reason why our 
data are consistent with the theoretical pole terms and why there is a difference in the value of v 2 . It should be 
noted that we are fitting over 7 orders of magnitude. An other independent source of discrepancy could be the tuning 
method and that the superfluous energy from the oscillon core goes out in shells 38] . These shells are observed to be 
present even for large time and their movement could increase the radiation power measured by our method. 

Let us turn our attention to the verification of the semi-empirical radiation law determined from the short-time 
evolution of various oscillon states. We shall confirm from a long-time simulation, that the evolution from an initial 
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(a)The theoretical and semi-empirical radiation law 



(b)i>2 as a function of e. 



FIG. 4: The measured data points, the theoretical radiation law and the results of the first fit with parameters to be found in 
Eq. (|61|l in the case of the (f> 4 theory. We see the points are approaching the theoretical curve as e decreases. 



oscillon state is driven by our radiation law. We claim two states to be identical if they posses the same e value with 
the same energy for a relatively long time. This should mean that the field configurations are very close to each other, 
the interaction with the radiation field is negligible. By testing these properties in two different simulations we can 
match these data and get a longer process. 

To compare the time dependence of the measured energy with the predicted one, we need the function E(e). This 
should not be a problem, as we can measure this function from the short simulations. The results of simulations are 
in very good agreement with the semi-empirical radiation law. We would like to emphasise that the radiation law is 
not a simple fit, it has been determined independently; we only use the initial e value and the function E[e) from the 
simulation to get the predicted curve. 

On Figure Owe see an e = 0.6 oscillon evolving. We used two simultaneously performed simulation for this graph 
(starting from e = 0.6 and e = 0.45), which perfectly fit together; another fact backing the hypothesis of adiabatic 
evolution. The solid lines follow the evolution observed in numerical simulations, while the dashed lines are the 
predictions of the semi-empirical radiation law (SERL). We end both curves at the same e value. 



2. Oscillons from kink-antikink initial data 



After being able to create and examine clean oscillon states we aim to identify, what the logarithmically decaying 
object is in our terminology that Geicke found in Refs. (35l. [39j. It turns out that the objects he found are composite 
oscillon states and there is a continuum of them, without any essential difference between these configurations in 
contradiction with what he conjectured. These objects do not obey the semi-empirical radiation law set for clean 
oscillon states in their early but very long period (« 150000) of life, they radiate more rapidly, as for the energy 
stored in the modulation modes of frequency has different radiation properties. However, their lifetime is in the same 
magnitude of oscillons with approximately the same frequency and energy. After their initial stage of evolution they 
begin to obey the semi- emp irical radiation law, their modulation, however, does not disappear. We shall confirm the 
logarithmic fits of Ref. [35j and in this stage of observation we can explain qualitatively what he found and how his 
results are related to the analytic considerations. 

Following Ref. [35j . we examine a kink-antikink pair initially at rest: 

+ 1 (66) 

d t <f>(x,t = Q) = 0. (67) 

For different values of a we examined the evolution of the initial state. The resolution was selected low (RES = 4) 
because we were not interested in the exact decay rate, we aimed to draw qualitative consequences. The results are 
collected in Fig. [7J For bigger values of a (a > 1.0) we find, that after an early period of formation approximately the 
same state arises. For moderate values of a we see different states arising in a relatively short period of time; these 
states radiate very slowly. Below a certain value of a (a ;$ 0.45) we do not observe any definite oscillon core. This 
behaviour is remarkable: we found the same pattern in the collision of a sine-Gordon soliton-antisoliton pair initially 
at rest, the same state arises for a > 1.4 and for the region a i; 0.8 we see no definite oscillon state. 



(x, t = 0) = tanh 



— tanh 
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FIG. 5: The evolution of an e — 0.6 and an e = 0.45 oscillon in the cj) 4 theory matched together; both the numerical data and 
the prediction of the SERL are plotted. The time evolution of the two oscillons are overlapping in a long time interval around 
t = 100000. 



To find out how these configurations radiate we ran a simulation with RES = 8, a = 0.8; this is the initial data 
of Geicke in Ref. [35| ■ The following figure shows us what the semi- empirical radiation law would predict and what 
really happens. The lumps emerging from kink-antikink initial data radiate faster then a clean oscillon due to its 
frequency modulation. The following asymptotic logarithmic fit considered in Ref. (35| works well: 

E(t) = = (68) 

W c + ln(t + V2- 10 5 ) V ; 

c= 0.52 ±0.12 B = 8.60 ± 0.07 . (69) 

Ref. (3f| uses m = \/~2, its result for the B parameter in the units we use in this paper is B g = 9.05, the reason for 
the difference between B and B g might be the use of Sommerfeld boundary condition in its simulation instead of 
compactification. We fitted the function which we get by extrapolating the theoretical result for infinitesimal, Eq. 
(|55|) oscillons to our case as well: 
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FIG. 6: The evolution of the kink-antikink initial data in the 4 theory for various a values; e versus time plotted. 



E(t) 



B 



B 
kit 

47T 



7.255 



theory — 

Bf it =8.190 ±0.007 



(70) 

(71) 
(72) 



We see that the numerical values from the fit and from theory are satisfactory close to each other. We note that such 
a logaritmical fit is not influenced by the value of in the theoretical radiation law, Eq. ([55)) . (Ref. [35[ claims that 
the results of its fit. B = 14.395 agrees with the theoretical result B = 14.503 in its units; however, the correct B 
value from theory with m = V2 mass scale is B = 10.260.) 
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FIG. 7: The evolution of the kink-antikink initial data in the 4> 4 theory with a = 0.8 parameter illustrated by the energy in the 
entire simulation. 



We conclude that the initial state observed by Ref. [35j evolves as a complex oscillon-like object. After the 'early' 
period of f ~ 150000 it looses a large part of its energy and an oscillon with modulated frequency is created. Through 
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its evolution the e value changes from 0.45 to 0.38. From the solutions of the semi-empirical radiation (SERL) law 
it can be clearly seen that the initial configuration decays much faster than an oscillon, however if we regard the 
t ~ 150000 as initial data, it obeys the radiation law with high precision despite the presence of modulation. (We 
did not plot the prediction of the SERL with to = 150000 in Figure [71 as the difference between the prediction and 
the numerical simulation results do not differ visibly.) It seems to us that the modulation degree of freedom is an 
adiabatically decaying mode. 

E. Verification of the theoretical radiation law 

We will examine oscillons in the specific symmetric <fi 6 theory in which oscillon radiation is the largest, i.e. when K 
is maximal. We determined the 175 value for this theory in subsection III El Thus, contrary to the </> 4 theory the value 
of K is maximal and the next to leading order corrections are significantly smaller compared to the leading order 
term, than in the case of the </> 4 theory because the (f> 6 potential is symmetric. These corrections originate in (j23|) . It 
should be kept in mind that because of the smaller pole term we cannot go as low in e as in the 4> 4 case. We collect 
the results of the simulations in Table I VIII in analogy to Table IVII 
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TABLE VII: Radiation power for oscillons with different e values in the symmetric 6 theory. 

We performed two fits: in the first we fitted both v 3 and b, in the second we used b = btheory an d fitted v% as we 
were interested in how accurately we could determine the value of ^3 from numerical simulations. 

(73) 

bfl = 0.89 ±0.1 (74) 

(75) 

btheory = 1 ■ (76) 

One can notice an anomalous point among the data points in Fig. [8l This anomalous configuration was further 
examined by us, but no visible anomaly in the field configuration was observed. The oscillon has a proper oscillating 
tail and a smooth core. Its radiation power was measured for different resolutions and also by using different values 
for the conformal factor k, but all these had no substantial effects. As it is, we do not know whether this anomaly is 
an awkward lattice effect or genuine one i.e. present in the continuum limit as well. 

We would like to investigate the problem of oscillon radiation from the perspective of field configuration as well. In 
Fig. [5]we show the field values as a function of time at x = 49.6 for the e — 0.352178 oscillon. We plot the theoretical 
prediction of this oscillating tail, so that the frequency of the wave is set to 3w, the amplitude is determined by theory 
and only the phase of the wave is fitted. We find satisfactory agreement in the case of wave amplitude and precise 
agreement in the case of frequency. If we compare the radiation power calculated from this plane wave and the one 
in Table IVTTl we find good agreement. Hence we conclude that the oscillons on the lattice loose energy via radiation 
and that out assumptions were correct when deriving the radiation law for small amplitude oscillons. 

There are various possible explanations for the discrepancy between theoretical and numerical results. Firstly, we 
could argue that the e values used in numerical simulations are too big and the theoretical calculation of the radiation 
amplitude only works for infinitesimal e values. For finite e values we can only expect an approximate agreement. 



dE 
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= -24\/2V3 exp 
= 0.4 ± 1.2 
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V3,theory = 0.9098 
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FIG. 8: The relation between the numerical data and the theoretical radiation law in the (f> 6 theory. We plotted the data points 
and the theoretically predicted curve only, as the results of the fits would only make the graph less understandable. 



Secondly, we cannot be sure about the initial data. Although the frequency of our objects is very stable there is no 
way to decide whether we work with undistortcd oscillons. The configuration with anomalous radiation emphasizes 
these problems. Lattice effects are less probable to play a role, as we see no major resolution and n dependence of 
the radiation powers. 
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FIG. 9: The field configuration observed in numerical simulations and the theoretically predicted wave amplitude for e = 
0.352178 oscillon in the 6 theory. 



We conclude that we can accurately determine the b parameter of the radiation law from the numerical simulation 
of oscillon decay in the case of both the (j> 4 and 6 case. The parameter v-i can be less accurately determined in the 
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case of the <fi 4 theory. There the next to leading corrections play a role because of the finite e oscillons and the value of 
v 2 is too small to be seen. In the case of the cj> 6 theory the next to leading order corrections are smaller and the value 
of ^3 is bigger, hence we can get the value of 1/3 from the simulations with satisfactory precision. The investigation of 
the radiation field of the oscillons show equally good agreement with the theoretical formulae giving a firm basis for 
the theoretical calculations from another perspective. 

IV. CONCLUSIONS 

In a general class of one dimensional scalar field theories we have computed the magnitude of the radiative tail of 
oscillons, determining their energy loss, in the small amplitude limit. The magnitude of the tail is non-perturbativcly 
small in the amplitude. We have used the Segur-Kruskal method of matched asymptotic expansion together with 
Borel summation techniques to calculate it. These results have also been verified numerically. We have also performed 
numerical simulations to compute directly the energy loss of oscillons, as well as the radiative tail. The numerical 
results are in a satisfactory agreement with the theoretical predictions. 
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